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to a finer lattice. We verify the existence of a pole mass as we take the continuum limit and deduce 
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I. INTRODUCTION 

The gluon propagator, although not an observable 
quantity, plays an important role in phenomenological 
non perturbative studies. A framework for such studies is 
provided, for instance, by the Dyson Schwinger equations 
(DSE) which recently have been applied, amoM other 
topics, to the study of the quark gluon-plasma ^ . This 
application is particularly important because it comple¬ 
ments experimental activity at RHIC and also because 
of the difficulties of applying lattice QCD to the case 
of non-zero chemical potential, although some progress 
is being achieved in this direction However in the 
Dyson-Schwinger approach the study of the gluon prop¬ 
agator is still inconclusive because of the various trun¬ 
cations needed to solve the coupled set of equations [Q. 
For instance, in the ghost-free axial gauge, many studies 
which used a simplified version of the three-gluon ver¬ 
tex supported an infrared enhanced gluon propagator of 
the form l/g^. Such a behaviour, driven by the vacuum 
polarisation diagram, was disputed by other DSE stud¬ 
ies 1^ as being due to a flaw in setting to zero the second 
scalar function that enters in the definition of the gluon 
propagator. Similar disputes also occur in the case of the 
Landau gauge where some studies, which assume dom¬ 
inance of the gluon-vacuum polarisation, find infrared 
enhancement Q whereas others [|j, which use a ratio¬ 
nal polynomial Ansatz for the self energies and vertices, 
find an infrared vanishing propagator. After including 
the ghost propagator, recent studies in the Landau gauge 
favour an infrared-vanishing gluon propagator [^. 

Lattice QCD provides a well defined theoretical frame¬ 
work for non perturbative physics and it is well suited for 
the study of the gluon propagator. A series of papers 0, 
appeared over the past couple of years, provide a de¬ 


tailed study of the behaviour of the gluon propagator in 
quenched lattice QCD in the Landau gauge. However, 
fixing to Landau gauge on the lattice is an iterative pro¬ 
cedure which stops upon reaching any of a large number 
of local minima, called “Gribov copies”. Since the global 
minimum itself cannot be reached, the effect of trading 
it for a local minimum is largely unknown. This unique¬ 
ness of the gauge condition puts us in a good position to 
obtain physical results on quantities like the pole mass 
which one may expect to be gauge invariant, and compare 
to the corresponding values used in phenomenology. The 
issue of gauge invariance of the gluon pole mass was ad¬ 
dressed in ref. For a certain class of covariant gauges, 
the Ward identities that determine the gauge dependence 
of QCD dispersion relations were derived. Using these re¬ 
lations one can examine the gauge independence of the 
poles of the gluon propagator. It was argued that if the 
loop expansion holds, then a non-zero pole of the trans¬ 
verse gluon propagator is gauge invariant. 

In ref. ||], we showed that a good description of the 
gluon propagator was provided by an Ansatz which ad¬ 
mits a dynamically generated gluon mass and thus 
points to an infrared regularised gluon propagator. By 
analytic continuation to negative values of we obtained 
an estimate of the pole mass. The existence of a gluon 
mass has important phenomenological implications [ pd] |. 
Total cross sections in hadron-hadron collisions, proton- 
proton elastic scattering and diffractive phenomena can 
be well understood if there is a finite correlation length 
for the gluon field . For instance in the Pomeron ex¬ 
change model of Landshoff and Nachtmann jl^ a gluon 
propagator which is infrared finite is shown to eliminate 
the troublesome singularity in the Pomeron calculation 
of hadron-hadron scattering. Whereas a bare gluon mass 
would lead to problems with unitarity, a dynamically gen¬ 
erated mass vanishing in the ultraviolet reproduces the 
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correct perturbative result for the gluon propagator and 
is consistent with unitarity. 

It is the purpose of the present work to check the ro¬ 
bustness of our earlier results on the gluon propagator as 
we take the continuum limit. We thus extend our previ¬ 
ous calculation on coarser lattices ||] to a finer lattice at 
/3 = 6.2. In ref. we included a comparison of results 
in different physical volumes demonstrating that, for the 
quantities of interest here, such as the pole mass and the 
zero momentum limit of the gluon propagator, a lattice 
size of about 1.5 fm suffices. Therefore for this study we 
use a lattice of spatial size ^1.7 fm. Our observable is 
the transverse part, D{q^), of the propagator. The excel¬ 
lent scaling behaviour which we observe for D{q^) enables 
us to extract accurately the change with /3 of the lattice 
spacing. In the infrared, it allows us to control the cutoff 
effects on the pole mass and to study its continuum limit. 

Our notation is the same as that of ref. ||^ and we refer 
the reader to [H for the details of our approach. 

II. SCALING 


Reasonable scaling was already observed in where 
we compared data at /3 = 5.8 and 6.0 on a lattice of size 
16^ X 32. Here we compare /3 = 6.0 and /3 = 6.2. 
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FIG. 1. Data at /3 = 6.0 (16^ x 32 lattice; crosses) and at 
/3 = 6.2 (24® X 48 lattice; triangles) fall on a universal curve. 

The results at /3 = 6.2 were obtained from 220 con¬ 
figurations generated by the UKQCD collaboration on a 
lattice of size 24® x 48. At /3 = 6.0 we used 200 coiifig- 
urations of size 16® x 32 from the NERSC archive [^ . 
Our Laplacian gauge condition consists of aligning along 
a fixed orientation the local 3-color frame built from the 
two lowest-lying eigenvectors of the covariant Laplacian. 
For implementation details, see 1^. Being now closer to 


the continuum limit, we find very good scaling behaviour 
for D{q^) as demonstrated in Fig. 1, where the two data 
sets fall on the same curve after applying the linear trans¬ 
formation 

\n{D/3=e,o{lnqa/3^e.o)) = ln(D/3=6.2(ln(ga^=6.2 - b)) -f c. 

The two fitted parameters take values b = 0.277 ± 0.022 
and c = —0.574 ± 0.053, which yields the scaling ratios 

= 1.02 ±0.14, = 0.758 ± 0.017. (1) 

This ratio of lattice spacings is consistent with that ob¬ 
tained from measurements of the string tension 
and very close to the value 0.729 obtained from the in¬ 
terpolation formula of the Alpha-collaboration for 
ro/a. Therefore, we use for the lattice spacing the value 
a“®(/3 = 6.2) = 2.718 GeV of Ref. |^], corresponding to 
a string tension ^/a = 440 MeV. 

The renormalised zero momentum propagator also ex¬ 
hibits good scaling. We obtain a value of (454(5)MeV)“^ 
in agreement with our previous value of ^ (445MeV)“^, 
both for a renormalization point of 1.943 GeV. 
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FIG. 2. The transverse gluon propagator multiplied by 
at /3 = 5.8, 6.0 and 6.2 in physical units. The fits to three mod¬ 
els are shown: Marenzoni (dashed dotted line) [17], Gornwall 
(solid line) [10], and model A (dashed line) of ref. [7] 

Our new transverse propagator data at [3 = 6.2 are 
shown in Fig. ^ together with our previous results. Since 
the /3 = 6.0 and 6.2 data agree so well over the whole 
momentum range, we can compare them to a variety of 
continuum models. We find that infrared enhancement 
like ( 9 ^)“^ 0, as well as Gribov-type infrared suppres¬ 
sion of the form D[q^) = Zq^l{q^ ± M^) L{q^,M‘^) |^] 
or Zq^/{q'^ + 2aq^-\-M'^)L{q^,M‘^) |^, are both excluded 
by our data. Here, the function L{q'^,M^) enforces the 
perturbative ultraviolet behaviour of the propagator: 
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L(g2,M2) = |iln[(q2 + M2)(q-2 + M-2)]| " . (2) 

The ansatz of Marenzoni et al. 

^(9 ) = ^(g2)l+a + (^2)1+0. ’ ( 3 ) 

with a non-perturbative anomalous dimension a, gives a 
better description of the lattice data (;^^/d.o.f = 18) but 
clearly underestimates the peak in Fig. A much better 
description is provided by model A of ref. 0 


D{q^) = Z 




_(g2 + M2)i+“ q‘^ + M‘^ 


L{q^,M^) 


( 4 ) 


with L{q^,M‘^) defined in (|^). Although this model 
gives the best fit to our data (x^/d.o.f ^ 2), it requires 
one more parameter, and lacks a theoretical motivation. 
Worse for us here, it cannot be analytically continued. 

Taking a different approach, Cornwall 0 allows for 
a dynamically generated gluon mass which vanishes at 
large momentum in accord with perturbation theory: 


D{q^) = Z 


{q^ + M^(q^)) In— 


4M2(g2)- 

"A2 


-1 


with 


M{q^) = M 


In [(q2 + 4M2)/A2] I 
In [4M2/A2] j 


( 5 ) 


As we already observed on coarser lattices, we find 
that also at /3 = 6.2 Cornwall’s Ansatz provides a rea¬ 
sonable fit to the data over the whole momentum range 
(x^/d.o.f = 7). We will thus keep it as one possible way 
to extrapolate to negative values of for the determi¬ 
nation of the pole mass. 


III. POLE MASS 

As explained in the Introduction, a phenomenologi¬ 
cally important question is whether the gluon propagator 
has a pole mass. The pole is a zero of the inverse propa¬ 
gator D~^{q^). We show this quantity in physical units 
on Fig. combining our data at /3 = 6.2 and /3 = 6.0. To 
determine the pole mass from D~^((f‘) = 0, an analytic 
continuation to negative q^ is needed. Given our finite 
amount of data, this extrapolation procedure is ambigu¬ 
ous. It is essential to test the robustness of a possible 
pole by comparing a variety of plausible extrapolations. 
Fortunately, the deviation of D~^{q^) from a linear func¬ 
tion of q^ in the infrared is small (much smaller than in 
Landau gauge |^), which increases confidence in the ex¬ 
trapolation. We consider here linear, quadratic and cubic 
polynomials in as well as Cornwall’s ansatz (i- All 
four choices indicate the presence of a pole. The variation 
in its location gives us crude systematic error estimates. 
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FIG. 3. Infrared inverse gluon propagator at d = 6.0 
(squares) and 6.2 (triangles). The various extrapolations to 
negative (f are: linear (dot-dashed), quadratic (dashed) and 
cubic (dotted) polynomials, and Cornwall’s model (solid line). 


Given the curvature of the data, a linear fit provides a 
lower bound for the pole mass, with a value of 505 MeV. 
This mass increases to 693(20) MeV with a quadratic fit, 
and reaches 800 MeV for a cubic fit. Note that the fitted 
coefficients of the higher powers of g^ keep decreasing, 
which indicates the soundness of an extrapolation based 
on a Taylor expansion. However our data are not suffi¬ 
cient to reliably fit higher-degree polynomials: in partic¬ 
ular, the results depend on the fitted momentum interval 
and on the presence or absence of data removed by the 
“cylindrical cut” of ref. 0 , whose purpose is to eliminate 
momenta most affected by lattice artifacts. In compari¬ 
son with these simple polynomial fits, Cornwall’s ansatz 
gives a pole mass of 669(6) MeV, consistent with that 
of the quadratic polynomial. The various extrapolations 
are shown together in Fig. 

By performing the same analysis at the three values 
of the lattice spacing a{P),/3 = 5.8, 6.0 and 6.2 we have 
studied, we can extrapolate the pole mass to the contin¬ 
uum limit. In Fig. ^ we compare such an extrapolation in 
a^, for pole masses obtained by fitting a quadratic poly¬ 
nomial or Cornwall’s ansatz. It can be seen that both give 
consistent results, at fixed lattice spacing as well as in the 
continuum. The continuum values are 632(38) MeV and 
592(14) MeV using the quadratic polynomial and Corn¬ 
wall’s Ansatz respectively. 
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FIG. 4. The pole mass at /3 = 5.8, 6.0 and 6.2 extracted 
from fitting to a quadratic polynomial (data shifted to the 
right for clarity) and to Cornwall’s Ansatz, as a function of 
the lattice spacing a^. The continuum values are obtained by 
linear extrapolation in a^. 


The reasonable robustness of our analysis, with respect 
to a change in the lattice spacing as well as in the an¬ 
alytic continuation procedure, leads us to conclude that 
there is strong evidence that a pole exists in the gluon 
propagator, with a mass of about 600 MeV. Because of 
the curvature of the inverse propagator, the systematic 
error in the extrapolation to negative is asymmetric. 
From the variation observed above with the choice of an¬ 
alytic continuation, we estimate it at ^ -30,-1-150 MeV. 
Based on our study of finite volume effects [|| , and given 
our present lattice size, we expect negligible finite size 
corrections to this pole mass. Comparing it to the glue- 
ball mass of 1.73(0.05)(0.08) GeV |19|, it appears close 
to one third of the glueball mass (rather than one half as 
sometimes speculated). A pole mass of 500 — 800 MeV is 
also within the range needed to fit experimental data in 
various phenomenological studies |11 

As a further, model-independent check on the value of 
the pole mass, we measured the correlator of the gluon 
field averaged over a time-slice, namely ® 


8 3 K K 

CW = (6) 

^ a=l ^=1 X X 


The exponential decay, if at all, of this correlator at large 
time is governed by the pole mass. Therefore, one can in 
principle measure the pole mass in a model-independent 
way. However, statistical noise limits the usefulness of 
this approach. We display the correlator in Fig.^. The 
data is of insufficient quality to measure an exponen¬ 
tial fall-off, and a fitting procedure is needed. Although 
this correlator is measured on the same configurations as 


D{q^), the various momenta are given a different weight, 
so that a fit to C{t) will give different results than a fit 
to especially after the cylindrical momentum 

cut in the latter. Since all momenta enter in C(t), we 
fit the only ansatz which describes the gluon propagator 
reasonably well over the whole momentum range, namely 
Cornwall’s model, directly to C{t) instead of D~^{q^) as 
before. The dashed line in Fig.gl shows the original fit 
of Cornwall’s ansatz to D~^{q^), which already provides 
a fair description of the data. The solid line represents 
a direct fit of the same 3-parameter ansatz to C(t), ex¬ 
cluding the first few time-slices which are contaminated 
by contributions from excited states. A simultaneous 
fit of the time-slice correlator data at (3 = 6.2 and 6.0 
yields a pole mass of 739(81) MeV, in agreement with 
the value 669(6) MeV extracted from D~^{q‘^), but with 
much larger statistical errors. 
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FIG. 5. Time-slice gluon correlator at /3 = 6.0 and 6.2. 
The dashed line shows Gornwall’s model fitted to 
after the cylindrical momentum cut; the solid line is a direct 
fit of the same model to the time-slice correlators, excluding 
the first few time-slices. 


If one attempts a model-independent determina¬ 
tion of the pole mass from the effective mass 
Weff(t) = —Ln{C{t + 1)/C{t)), one obtains a value of 
702(163) MeV, poorly determined but consistent with 
the direct fit. Therefore, our different analyses give pole 
masses ranging from 592 MeV (continuum extrapolation 
of pole of propagator fitted by Cornwall’s Ansatz) to 
739 MeV (fit of time-slice correlator). Taking into ac¬ 
count the asymmetry of potential systematic errors, we 
estimate the gluon pole mass to be ^ 600^30° MeV. 
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IV. CONCLUSIONS 


[81 C. Alexand rou, Ph. de Forcrand and E. Follana, |iep- 


lat/0009003; Phys. Rev. D 63 (2001) 094504. 


We have extended a previous lattice study of the gluon 
propagator in the Laplacian gauge to a finer lattice and 
found good scaling behaviour. We confirm the existence 
of a pole as we approach the continuum limit. Applying 
a variety of different hts we extract a pole mass in the 
range of 600 — 30 + 150 MeV in accord with the value 
found in phenomenological studies for the description of 
hadron-hadron scattering. 

It would be very interesting to substantiate the gauge- 
invariance of the gluon propagator pole mass by simi¬ 
lar studies in other gauges. Ref. |^| proposes a non¬ 
local, gauge-invariant gluon propagator based on the 
(Coulomb-like) Laplacian gauge in three dimensions. It is 
argued there that the pole mass of this propagator deter¬ 
mines the vector- pseudoscalar- mass splitting My — Ms 
in heavy quarkonia. This relation is consistent with the 
SU{2) numerical results presented. For QCD, using the 
experimental values for the splitting in cc and bb systems, 
the implication is that the pole mass is ~ 420 MeV. This 
is somewhat low compared to our estimate. In Landau 
gauge, our first attempt ||] showed no indication for a 
gluon propagator pole. However, lattice Gribov copies 
might play a crucial role there. Moreover, Landau gauge 
appears more sensitive than Laplacian gauge to finite- 
size effects. Thus this question requires larger lattices as 
studied in |^. 


[9] R. Kobes, G. Kunstatter and A. Rebhan, Phys. Rev. 
Lett. 64 (1990) 2992; Nucl. Phys. B355 (1991) 1. 

[10] J. M. Cornwall, Phys. Rev. D 26 (1982) 1453. 

[11] G. Paris! and R. Petronzio, Phys. Lett. B 94 (1980) 51. 

[12] F. Halzen, G. Klein and A. A. Natale, Phys. Rev. D 47 
(1993) 295; M. B. Gay Ducati, F. Halzen and A. A. Na¬ 
tale, Phys. Rev. D 48 (1993) 2324. 

[13] P. V. Landshoff and O. Nachtmann, Z. Phys. G 35 (1987) 
405. 


[14] http://qcd.nersc.gov 


[15] G. S. Bali and K. Schilling, Phys. Rev. D 47 (1993) 661. 

[16] M. Guagnelli, R. Sommer and H. Wittig [ALPHA collab¬ 
oration], Nucl. Phys. B 535 (1998) 389. 

[17] V. N. Gribov, Nucl. Phys. B139 (1978) 1. 

[18] P. Marenzoni, G. Martinelli, N. Stella and M. Testa, 
Phys. Lett. B 318 (1993) 511; P. Marenzoni, G. Mar¬ 
tinelli and N. Stella, Nucl. Phys. B455 (1995) 339. 

[19] C. Morningstar and M. Peardon, Phys. Rev. D 60 (1999) 
034509. 

[20] O. Philipse n, Phys. Lett. B521 (2001) 273 ; hep- 


lat/011011+ 


Acknowledgements: The 16^ x 32 lattice configurations 
came from the Gauge Gonnection archive ||T^, and the 
24^ X 48 were provided by the UKQGD collaboration. 


[1] G. D. Roberts and S. M. Schmidt, 
Phys. 45S1 (2000) 1. 

[2] Z. Fodor and S. 


Prog, of Part. Nucl. 


D. 


hep-lat/0106002|. |hep-lat/0110102 

[1 


Katz, 


ration, |hep-lat/0107002| , |hep-lat7C)110223| 


QCD-T ARO collabo- 
S. Ejiri et ah. 


hep-lat/011008C. 


[3] G. B. West, Phys. Rev. D 27 (1983) 1878. 

[4] S. Mandelstam, Phys. Rev. D 20 (1979) 3223; N. Brown 
and M. R. Pennington, Phys. Rev.D 38 (1988) 2266; 
Phys. Rev. D 39 (1989) 2723. 

[5] U. H’abel, R. Konning, H. G. Reusch, M. Stingl and 
S. Wigard, Z.Phys A 336 (1990) 423;ibid 435; M. Stingl, 
Z. Phys. A 353 (1996) 423. 

[6] L. von Smekal, A. Hauck and R. Alkofer, Phys. Rev. Lett. 
79 (1997) 3591; Annals Phys. 267 (1998) 1. 

[7] D. Leinweber, J. I. Skullerud, A. G. Williams, and 
C. Parrinello, Phys. Rev. D 60 (1999) 094507; Erratum- 
ibid. D 61 (2000) 079901; F. D. R. Bonnet, P. O. Bow¬ 
man, D. B. Leinweber and A. G. Williams, Phys.Rev. D 
62 (2000) 051501. 


5 



























